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1 Abstract 

Self-consistent inouvement of initial perturbation in density, velocity and grav- 
itation potentail on the background of the stationary cylindrical configuration 
of the gas with gravitation and pressure in Lagrange variables have been stud- 
ied. The nonlinear partial differential equation for description radius motion has 
been obtained. The linearization of this equation is reduced to a Klain-Gordon 
equation which has an analytical solution. 

2 The Model 

Self-consistent mouvement of initial perturbation in density, velocity and gravi- 
tation potentail on the background of the stationary cylindrical configuration of 
the gas with gravitation and pressure in Lagrange variables have been studied. 
In [0 it is shown that a cylindrical symmetry admits an equilibrium state of the 
gas with selfgravitaiting and pressure. The cylinder is supposed infinity long 
along the axes and all physical values are dependent on the radius of cylinder 
only. The equilibrium is provided by the equality of gravitational force and force 
of gas pressure in every point. In this article the model used in [1] is appeared as 
a background for a small perturbations in velocity, density, pressure and grav- 
itataitional potentail. In general case the perturbation may depend from two 
coordinats and have very complicated geometrical form. It is not possible to 
study two-dimensional nonlinear motion in general case. But the small devia- 
tion from the stationary state admits the principle of superposition. According 
to this it is possible to study radius and longitudinal motions independently 
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and then to sum them taking into considearation the fact that they are vectors. 
This article is dedicated to study pure radius motion. 

The model of perturbation depending on radius of cylinder only is choosen. 
It means that the cylinder infinite along the axes of symmetry makes (as a 
whole) radius motion by the influence of initial perturbation of stationary state. 
The perturbation is distributed along the axes of cylinder thus its value is de- 
pendent on the radius coordinat only. Two causes are able to generate the initial 
perturbation: 1) the velocity is not equal to zero for particles in the equlibrium 
state or 2) displacement of a particle from the equilibrium point is not equal to 
zero. The combination of the both is possible as well. 

The mathematical description of this model is represented by the Cauchy 
problem for three-dimensional nonstationary partial differential equations of 
motion, continuity, Poisson equation and algebraic equation of state: 

dv VP 

^ + {v.W)v^-—-W<f (1) 
ot p 

|f + V-M = (2) 

= iirGp (3) 

P = Ap'' (4) 

whis initial conditions 

v\t=o = ^(^^ 0), P|,^o = P(f, 0), $1,^0 - 0), = p(f, 0) (5) 

where P is pressure, v - speed, $ - gravitational potential, p - gas density, t - 
time, A = const, 1 < 7 < 2. v{f,0), P(r, 0), $(r, 0) and p(r, 0) are specified 
functions. The system (|^) - describes hydrodynamic motion in the ideal 
classic gas with selfgravitation and pressure. 



3 Characteristic Values and Parameters 

Characteristic values are needed to transfer from dimension to dimensionless 
equations. A choise of characteristic values is defindent on the physical model. 
Having been choosen these values constitute dimensionless coefficients of the 
equations. Depending on the initial and/or boundary conditions these coeffi- 
cients are able to play a role of small or big parameters. There are two kindes 
of characteristic dimensional values in this problem. One kind of them is con- 
nected with stationary gas configuration (the background) and the second one - 
with perturbations expanded on the stationary background. Some various per- 
turbations dependent on the initial conditions have distinction by different time 
characteristics (a period, for example). To compare them it is essentially to 
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have a time scale not dependent on the perturbation. There is only one choice 
to do so - to attach the time scale to the background. 

The stationary configuration of the gas will be discribed by a series of the 
characteristic values - Jeans lenght Lq, gravitaitional potential $0, gas pressure 
Po, mass of particle m, concentration of particles A^o, density po, themperature 
Tq, sound speed cq: 

Ll = ^ <i>o = ^ Po^iVofcTo c^ = A^p,-^-^ (6) 
that is why: 

t scale — 1 ^ C^) 

Full enegy of the gas consists of gravitaitional energy and energy of the gas 
thermal expansion. The full energy is distribited between these two components 
depending on the adiabatic index 7. The function /i(7) dicribes this distribution 
as follows: 

.(7) - (8) 

Equations (|^) - (||) give the correlation between a number of particals in the 
cube with rib length equal to Jeans length and parameter /i(7): 

NoLl = ^7 /.(7). (9) 

The perturbation will be characterized by the following values definded by the 
initial conditions: the length of wave Ao, period of perturbation to, velosity vo\ 
characteristic Mach number and characteristic dimensionless wave length. They 
are: 

Ao = voto Mq = — Ko = (10) 

Co ^0 

A parameter a will also be used: 

to 



From (|), and (|Tl|) it follows: 



(11) 



(12) 



4 The Initial System of Equations in the La- 
grangian Variables 

To convert from dimesionall equation (|l]) - to dimesionless ones the di- 
mensionless radius, density, pressure, gravitaitional potentail and velosity are 
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introduced according to the following rules: 

£, = — S=— P=-^ '^=;sr v==— . (13) 
Lq Pa Pq $0 vo 

In the cylindrical system of coordinates the dimesionless radius component of 
equations (|^) - (^) are 

dv dv 02 dP dip 

a7 + «iV3^ = -Ta^-«3^ (14) 

''U^)-a,5 (16) 



P = ,5'' (17) 

with initial conditions 

v|i.o=v(e,0) p|,^o=P(e,0) <^l4.o = <^(e,0) <5|i.o = ^(e.O) (18) 

where v(^,0), p(^,0), (/)(^,0) and (5(^,0) are specified functions and 

voto Poh $0^0 inGpoLl 

ai = — — 02 = 03 = 04 = . (19) 

Lq Po^oVq ^QVq 9o 

According to the (||) and the expressions for oi - 04 become: 

aMo aMo a . 2 / n ,on\ 

oi = — ^ 02 = — p 03 = ^ — rnrr "4 = 47r 7Ai(7). (20) 

Finally in Euler variables the equations being looked for are obtained by ex- 
cluding the a from the (|20|) and substituting the result into the (llj) - (|l6|): 



1 5v 9v 1 dP 1 



7|:f^t)=4^V(7)^ (23) 



with initial conditions ( |l8| ) . But the problem under consideration may be solved 
in the Lagrange coordinates only. Applying the transformation rules ( (|] ) : 
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(cr is a dimensionless Lagrange mass variable) to eqations - gives the 
system of equations describing the problem under consideration in Lagrangian 
system of variables: 



1 dv dp 5^ dip 

kq Ot 7 da jfi{j)MQ da 



with initial conditions: 



(25) 



— + Ko<5' 4-^=0 (26) 
OT da 

9 f,2,d^\ _ .2, 



<=o^H^^0) P\t=o^Pi^^O) 0|,=o = '^(^'O) S\t^,^S{a,0) 

Clt=o = e(a,0) 5(0,0) = 1 (28) 

where v(cr, 0), p{a,0), (f>{a,0), S{a,0) and ^(cr, 0) are specified functions. The 
dependence of function ^(cr, 0) from a will not be used in this article and, that 
is why, will not be studied. 



5 The Development of the Radius Motion Equa- 
tion 

In this section the system of equations (^5|) - (|2^) will be reduced to one equation 
for Euler coordinate ^(cr, r). The dependence on (cr, t) will be omitted expect 
the case r = 0. It means that d means S{a, r), but for r = we will write S{a, 0). 
The Poison equation ( p7| ) may be integrated directly with initial condition: 



dip 
d^ 



= const. (29) 

cr=0 



Because of the equation (|27| ) structure const should be finite but it is not es- 
sential and it doesn't matter which the const is. The first integral of Poisson 
equation is: 

= 47r7Ai(7)a - L(t) (30) 

da 

where L{t) is some function will be defined later from the conditions of inde- 

a<T 



pendence all physical values from the axes of symmetry. Excluding the S^S 



from (|27|) and (gg|) the equation: 

g 9v ^ e dP 4n^jfiij)a ^ Ljr) 
kq dr 7 da 7^(7)Mq 



is obtained. 
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The foUowing transformation is connected with the equation of continuity 
6|). Using the definition of velosity in the Lagrangian variables: 



and dividing (26) by 5 the equation: 

d 1 1 92^2 



dr 5 kq drda 

is obtained. Its integral is: 



(33) 



-6 ~ Y 97 = + ^ ~ Y 97 ^^^^ 



The function G(cr) will be defined from the initial conditions. From the defini- 
tion of a (p3) and equality: 



it follows that: 



Due to this 



da '^=0 da 
da '^=0 5{a,0)' 



(5(cr,0) 2 9a ~ (S(cr,0) 
Finally, the integral of continuity equation is: 

I Ko 1 — Ko 



(35) 
(36) 
(37) 



S{a,0) 2 9f7 d{a,0) 

Let's designate 



G(») (38) 



The next step is substitution the p from ( p8| ) into the equation of state ( p7[ ) 
and motion equation (25): 

C 92^ 5 fnode 4^V(7)^-L(r) 



Ko dr^ 7 9(7 \ 2 9(7 y jp,{j)Mq 

The function L(r) is found out from the follows condition: all physical values 
are not dependent on the axes of sylinder. Let's use it for a = ^ = 0. It is foUws 
from (g^): 

L(r) - (41) 
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The function G{a) is found out from the initial conditions (|2^). To calculate 
it we introduce a parameter ^o(o'): 

^ -CoM (42) 



5t2 



(<7,0) 



Using - (|39|) and calculating (^) for r = the equation for function G{a) 
is obtained: 



G{a) ~ - G{a) + - - ---pr - 77 = where ' = -. 

(43) 

To find out the initial condition for this equation we calculate the ( p8| ) at point 
(7 = and obtain: 

G(0) = 0. (44) 
So, the equation of radius motion which should be obtained in this section is: 

e d [node , ^^^^^^ 

with initial conditions 



= e(a,0) (46) 

t=o 



where v((t, 0), p(cr, 0), (/)(cr, 0), (5(o', 0) and ^(cr, 0) are specified functions. Func- 
tion G{a) is definded by the equation ( p3| ) with initial condition (^). The 
equation ( p6| ) is nonlinear nonstationary partial differential equation which can- 
not be solved analytically in general but only after linierization. 



6 The Linearization of Radiuse Motion Equa- 
tion 

The sections 5 and 6 ar edevoted to studying the small perturbations of the 
gas. An arbitrary radius motions are described by the equation ( ^ ) with initial 
conditions (|46|). According to the Model the transformation from (^) - ( ^6| ) 
to linear approximation will be made now. Let's assume that a particle is in 
an equilibrium in the point of radius with Lagrange coordinat ctq- Because 
of smallness of deviation the particle from the point of equlibrium the Euler 
coordinat of the eparticle may be presented as a sum: 

^{a,T)^^a + ^{<^.r) where |V'('^, ^)l < and Co = ^(^^o, 0) (47) 

To study a generaly case let assume that the particle is displaced at the moment 
of time T = from the poin of equilibrium to new poin with lagrange coordinat 
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a — <7o + Act and has a velosity equal to v(f7o + Ac). So, the initial conditions 
according this model are 



^(a,0) =e('To,0)+V(^,0) 



0(a,O) 



(48) 



where 



l/'(cr,0) = *(5[cr- ((To + Act)] 1^(cr, 0) = *(5[cr - ((Tq + Act)] (49) 



where 4' and ^' are parameters of the problem, and d[a — (ctq + Act)] is a delta- 
fuction. Let's designate: 

CT = CT + Act (50) 

In addition to (^^ the condition 



node 
2 da 



< S(ct) 



(51) 



will be used for linearisation of the equation (^5|). (^Tj) allows to simplify 



E(ct) 



«S(a)^+^ + (7 + l)y^S^(^) (52) 



2 Sct 

and after substitution the (|5^) into (|45|) the foUwing equation is obained: 



2 c)ct2 



(53) 



where S' = To linearized ( |53| ) use the (^. This substitution gives the 
equaton 

a^V' KqCo ^^V- 47r2K29^ + 47r2KoCT kq^oK 



EoAf2 da Mi^o 



-^7+1 



= 



To simplify this equation denote 

„,,2 _ '>'0?0 



V-7 + 1 
^0 



_ 47r2(7 + l)«2^ So 



SoM2 



^0 



47r2KoCT 



Then the equation (54) become as follows: 



(54) 



(55) 



(56) 



with the initial conditions ( p9[ ) . i? has a sense of acceleration. According to this 
two new characteristic values are introduced: 



a* = ^ and r* 
B 



CT" 

Wo 



(57) 
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The equation (pq) will be transformated now by the substitution 



V'(ct, t) = u{a, t) exp + C (r*)^ 

where u{a, r) is a new function. The result of the substitution is: 







(58) 



(59) 



For this equation the Cauche problem is calculated with the following initial 
conditions which are obtained by transformation according to (pq): begin: 



u{a, 0) = V'(cr, 0) 



du{a, T 



ijj{a, 0) 



(60) 



T=0 



It is the Klaine- Gordon equation ([g|) 



7 The Solution of the Klaine-Gordon Equation 

To solve the Kline-Gordon equation the Riemann method of the Cauchy problem 
solution will be used (jj]). The equation ( ^ ) is a hyperbolic equation. In the 
canonic coordinates (i^, 77): 



2a*wo 



V = 



a — WqT 
2a*wa 



the equation (59) is reduced to the form: 







with initial conditions ( |60|) transformated according to the rules: (|6l]): 

du{a, t) 



(61) 



(62) 





/ du 




2a* 




di] J 



y—il 



dr 



T = 

(63) 

In the Riemann method the solution of the equation ( |6^ ) is determined by the 
(|63[) and Riemann function V{v, rf) satisfying the conjugatee equation: 



dvdrj 4 



(64) 



with initial conditions: 



(65) 



9 



Let's find out the Riemann function V in the form 
V = N{p) 



where p — [v — i')[f} — if) (66) 

wherer [v^ fj.) is the coordinate of the point where the solution found out and 
{vji}) belong to the line in which the initial conditions are given. From (67) it 
follows that r = correspond to ;> = 77. The substitution ( |67| ) into the (3^ 
gives the ordinary differential equation: 



n" + -N' + N{p) = 
P 



which has a finite solution: 

N{p) = Jo {V{^-i>)ifl-v)) 



(67) 



(68) 



where Jq is a Bessel function of null order. The initial conditions (^5|) are 
satisfied. Let's return to the equation (|62|). According to the Riemann method 
the solution of the equation (^2|) satisfying the conditions ( |63| ) has a form: 

dO — 



/ du 









V ~ rj 



y Jo {\/ {v - - rf)^ 



dv 



(69) 



Let's transform now the initial conditions 
dependence on v: 



\j from the dependence on a to 

a*wo{v + = u{a, 0) ^-^5 [a- (7(0)) = 









^5 


a 


*Wo(i> 


+ fl)\0=f) 


-a(0) 


/ du 


du 
drj 




du{a, t) 
dr 


T = 


*(5(cr- 


am = 










a 


*Wo{l' 


+ fl)\o=f, 





1 

2a* 



After substitution the ( |70| ) and 

T 



v-ri= — [v - v){i> -rj) = {wot) - (ct - cr(0) ) 

a 



(70) 



(71) 



into (|69|) the solution of the equation ( |62| ) with initial conditions (|60| ) obtain: 
u(cr, r) = i {(5 [woT +[a -a)]+5 [wqt - [a - a)]} + 
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' I r ( VK,r)^-(^-<T(0))^ ^ 



2wo 

This formula together with (|58D define the solution of linear problem 



8 Results and Discation 

The obtained solution has a formal character of oscillations. But there are two 
reasons due to which this conclusion should be triated critically. First: This is 
a solution of a linear equation. The area of the initial conditions where linear 
solution differs slightly from the nonlinear one is not discuss in this article. The 
second: The solution has been obtained in the Lagrangian variables but when 
we are speaking about oscillation Euler variables are meant. The transformation 
from Lagrange variables to Euler sones is not performed in this article also. 
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